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The excitation of the degenerate Ei carbonyl stretching vibrations in dimanganese decacarbonyl is shown to 
trigger wave packet circulation in the subspace of these two modes. On the time scale of about 5 picoseconds 
intramolecular anharmonic couplings do not cause appreciable disturbance, even under conditions where the 
two El modes are excited by up to about two vibrational quanta each. The compactness of the circulating 
wave packet is shown to depend strongly on the excitation conditions such as pulse duration and field strength. 
Numerical results for the solution of the seven-dimensional vibrational Schrodinger equation are obtained for 
a density functional theory based potential energy surface and using the multi-configuration time-dependent 
Hartree method. 



I. INTRODUCTION 

Recently, considerable attention has been paid to the 
excitation of degenerate electronic and vibrational states 
using circularly polarized laser light. Manz and cowork- 
ers have been the first to focus on the possibility of the 
excitation of electronic ring currents by laser-triggered 
formation of time-dependent hybrid states, e.g. in a Mg- 
porphyrin modeP^ as well as in the diatomic AlCl.'^ A 
particularly interesting aspect has been the accompany- 
ing generation of giant magnetic fields."^ Subsequently, 
this work was extended to the vibrational domain in the 
electronic ground state. Here, degenerate bending vibra- 
tions of linear triatomic species such as -'^^^CdH2^ and 
FHF~ have been excited to yield unidirectional pseu- 
dorotations. In contrast to previous work on the pseu- 
dorotation of Nas in the electronic B state excited by 
linearly polarized light, ^ unidirectionality is achieved by 
circular polarization. 

The present contribution is motivated by the ques- 
tion whether the dynamics of degenerate vibrations 
can be driven in more complex molecules, where one 
has to compete with intramolecular vibrational energy 
redistribution.!^ Promising candidates are metal-carbonyl 
compounds due to the strong and spectrally distinct in- 
frared (IR) absorption of the carbonyl stretching vibra- 
tions. Previous studies of metal-carbonyl compounds 
have been focused especially on laser-driven CO vibra- 
tional ladder climbing .'^'^1^' A particularly interesting ex- 
ample has been the experimental demonstration of lad- 
der climbing up to the vqq = 6 state in carboxyhe- 
moglobin by JofFre and coworkers .'^i' For this example, 
Meier and Heitz have developed a non-reactive gas phase 
model to simulate the excitation dynamics by using local 
control theory.E^ It should be noted, however, that the 
main goal of the experiment had been to achieve con- 



trolled metal-carbonyl bond breaking assisted by anhar- 
monic couplings between the initially excited CO mode 
and the metal-CO bond. A quantum dynamics study 
employing a reactive model provided support for the ex- 
perimental finding that this anharmonic coupling is too 
weak to ca use bo nd breaking on the time scale of a few 
picoseconds .Em The lack of sufficient anharmonic metal- 
CO coupling has also been found in the related model 
study of MnBr(C0)5 using a three-dimensional reactive 
Hamiltonian.ES Further experimental progress in vibra- 
tional ladder climbing has been reported in Ref.'16'where 
an evolutionary algorithm was used for pulse shaping 
such as to highly excite a CO mode in W(C0)6. 

In the following we will consider dimanganese de- 
cacar bonyl Mn2(CO)io (cf. Fig. [l]) as a model sys- 
tem. In the past it was in particular the photochem- 
istry of Mn2(CO)io that attracted most attention (see, 
e.g., Ref. |T7] and references cited therein). More re- 
cently, the prospect of having ten coupled carbonyl vi- 
brations attracted the interest of two-dimensional IR 
spectroscopists.'^SHMJ This included in particular a de- 
tailed investigation of the coupled CO fundamental and 
low-order overtone and combination tone vibrations us- 
ing anharmonic vibrational perturbation theory for the 
subset of CO modes, which were treated by anharmonic 
couplings up to the fourth order.l^The anharmonic level 
shifts for all CO vibrations were found to be below 10 
cm^^ . 
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In Section |II A| we will present results on the quan- 
tum chemical characterization of the global minimum 
structure of the Mn2(CO)io model system. Based on 
this structure we will develop a seven-dimensional (7D) 
model Hamiltonian in Section |II B[ which is designed 
such as to describe the dynamics of the two degenerate 
El CO vibrational normal modes and their anharmonic 
coupling to the remaining degrees of freedom, which will 
be called a bath. Further, we will comment on the nu- 
merical solution of the time-dependent Schrodinger equa- 
tion, which is done using the multiconfiguration time- 
dependent Hartree (MCTDH) approach to wave packet 
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FIG. 1. Optimized D4d structure of Mn2(CO)io using 
DFT/BP86 with a TZVP basis set. 



dynamics. I^^l^ Numerical results on the laser-driven dy- 
namics will be discussed in Section IIIII and conclusions 
are present in Section [IV[ 



II. THEORY 

A. Characterization of the Equilibrium Geometry 

The structure of Mn2(CO)io (cf. Fig. [I]) has been op- 
timized under symmetry constraint employing Den- 
sity Functional Theory (DFT) with the BP86 functional 
and a triple zeta (TZVP) basis set as implemented in the 
Gaussian 03 program package.'^S Available crystal struc- 
ture data are well reproduced as seen from the compari- 
son given in Tab. I. 

The stationary points were further characterized by a 
normal mode frequency analysis. The results for the IR 
active CO stretching vibrations are summarized in Tab. 
II. Overall, the agreement with experiment is quite rea- 





calc. Ref . Ref . \M Ref. ITfl 


Mn-Mn 


2.97 2.92 2.98 2.90 


Mn-COax 


1.81 1.73 1.80 1.82 


Mn-COoq 


1.85 1.83 1.87 1.86 


C-Oax 


1.16 1.16 1.15 1.15 


C-Oeq 


1.15 1.16 1.15 1.14 


Mn-Mn-COcq 


86.8 86.2 86.6 86.1 



TABLE I. Comparison between calculated geometrical pa- 
rameters of the D4d stationary point of Mn2(CO)io and crys- 
tal structure data (bond lengths in A, bond angles in degrees). 



sonable even at the level of the harmonic approximation. 
Moreover, comparison with Ref. [T9l indicates, that there 
is only a small dependence on the basis set (double vs. 
triple zeta). 

B. Model Hamiltonian and Equations of Motion 

The potential energy surface has been expressed in 
terms of normal mode coordinates, i.e. V = V{Q,q). 
Thereby, the two normal mode displacement coordinates 
describing the degenerate Ei modes Q = (Qi, (32), which 
are treated beyond harmonic approximation, have been 
separated from the remaining 3iV — 8 harmonic bath 
modes q. The coupling, V'c(Q,q), between both sets of 
coordinates is modeled in linear order with respect to q. 
Neglecting any contributions from global rotations the 
molecular Hamiltonian becomes 

i?mol = i?Q+^fq + K(Q,q), (1) 

with the part describing the anharmonic coordinates 

HQ = ^{P? + Pi) + V{Q,q = 0), (2) 
the harmonic bath part 

(3) 

i 

and the interaction 

V-c(Q,q) = '^^Q,^^"^ q--F(Q)q. (4) 

For Mn2(CO)io there are 58 bath modes, but not all of 
them will couple appreciably to the Q modes. Using the 
forces F(Q) the reorganization energy has been calcu- 
lated according to the assumption that the frequencies 
do not change, i.e. E^"""'^^ = [F,(Q)/a;,]2.ISl This quan- 
tity has been averaged on the Q grid and the result is 
used to select the most relevant bath modes. The cut-off 
for the selection has been chosen such that modes whose 
averaged reorganization energy is substantially smaller 
than their harmonic frequencies are neglected. This gave 
a total of five relevant bath modes, which are combined 



mode 


polariz. 


calc. exp.l2SI exp.l22l calc.li2l 


B2 


z 


1988 (788) 1983 1992 1981 


El 


(x,y) 


2014 (2224) 2014 2025 2005 


B2 


z 


2045 (1348) 2045 2053 2036 



TABLE IL Comparison between calculated (harmonic, not 
scaled) and measured (Ref. [22] in gas phase, Ref. [SS] in n- 
hexane) frequencies of IR active carbonyl stretching vibra- 
tions (in cm~^ ). Calculated IR intensities are given in paren- 
thesis (in km/mol). 
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TABLE III. Most important bath modes according the aver- 
aged reorganization energy (average over the full Q 
grid and - in parentheses - over a reduced grid from -1 to 1 
ao(a.m.u.)^^^). The parameters of the MCTDH implementa- 
tion (A^dvr: number of DVR points, A^spf: number of single 
particle functions) are given as well. Notice that modes q2 
and gs as well as 53 and 54 have been combined. 



mode 


(cm-1) 


(cm-1) 


grid 
(ao(a.m.u.)i/2) 


A''dvr 


A'"spp 


qi (Ai) 


425 


468 (72) 


-7.5:7.5 


45 


10 


q2 (Ai) 


1999 


822 (122) 


-1.7:1.7 


25 


10 


<Z3 (E2) 


2016 


1426 (204) 


-1.7:1.7 


25 


10 


<?4 (E2) 


2016 


1649 (239) 


-1.7:1.7 


25 


10 


95 (Ai) 


2106 


4646 (692) 


-1.7:1.7 


25 


10 



with the anharmonic coordinates to a seven-dimensional 
model. The bath mode parameters are compiled in Tab. 



Ill 



Notice that the average reorganization energy de- 
pends, of course, on the grid size. In the table we report 
the average over the full grid as well as over a reduced 
grid onto which the wave packet will be mostly localized. 
The major effect on the dynamics of the Q modes will 
come from the coupling to the other CO vibrations, not 
only because the forces are the largest, but also because 
of the near resonance. Among the lower-frequency modes 
only qi couples appreciably; it is a breathing mode of all 
equatorial CO groups with respect to the two Mn centers. 

The molecular Hamiltonian is supplemented by the 
molecule-field interaction that is described in dipole ap- 
proximation 



-fffieidW = -d(Q)E(i), 



(5) 



where the dipole moment surfaces for the different po- 
larization directions have been obtained along with the 
potential energy surface for the selected anharmonic de- 
grees of freedom. Notice that within the approximation 
given by Eq. Q there is no coupling to the other two 
IR active carbonyl modes (cf. Tab. II). Moreover, the 
latter are polarized in z-direction and will not be directly 
excited by the circularly polarized laser pulse. For the dy- 
namics simulations it will be assumed that the molecule 
has been pre-oriented with the coordinate system as in- 
dicated in Fig. [1] The laser field will be taken to be of 
the form 




= EoQ{t)Q{tp - t) sin2(i7r/ip) 




(6) 



will be chosen such as to achieve the excitation of a vi- 
brational superposition state with respect to the two Ei 
modes, which corresponds to a wave packet performing 
a clockwise circulation on the potential energy surface 
F(Q, q = 0). The respective normal mode displacements 
are shown in Fig. [2] Close inspection of this figure shows 
that during this circulation, the antisymmetric combina- 
tion of stretching motions of equatorial CO groups which 
are opposite to each other " moves" around the symmetry 
(z) axis. Of course, this classical picture applies only as 
long as the wave packet stays compact and an important 
question to be addressed will be the regime of validity of 
this classical analogue. 

The 7D time-dependent Schrodinger equation 



[ifmoi + i?ficid(t)]*(Q,q;t) (7) 



Here tp is the total pulse duration, uj the carrier fre- 
quency, and Eq the field strength. These parameters 



has been solved using the MCTDH approaclPmi 

as im- 
plemented in the Heidelberg program package.'^ The 
seven-dimensional wave function ^'(Q, q; t) is represented 
on a grid in terms of a harmonic oscillator discrete vari- 
able representation (DVR). A total of 25 DVR points on 
the grid [—1.7 : 1.7]ao(a.m.u.)^/^ together with 12 sin- 
gle particle functions (SPFs) are used for each of the Q 
modes which are treated as a combined MCTDH parti- 
cle. The parameters for the bath modes are compiled 
in Tab. |III[ Potential energy and dipole moment sur- 
faces have been fitted on this grid using the POTFIT 
algorithm .1221 Using this set of SPFs the natural orbital 
populations have been below 10~^. The actual propaga- 
tion is performed using the variable mean field scheme 
with a 6th order Adams-Bashforth-Moulton integrator. 
The vibrational ground state, |^'o), has been calculated 
using the improved relaxation schemd^ some uncoupled 
(one-dimensional zero-order) states in the Q-subspace 
will be used for the discussion of the dynamics. 



III. RESULTS AND DISCUSSION 

In the following we will present results for the time- 
dependent wave packet dynamics which have been ob- 
tained for a carrier frequency w, fixed to the funda- 
mental transition of the Ei modes, which is at 2014 
cm^^ . The total propagation time has been set to 
5 ps. First, we will discuss the effect of changing 
the field strength Eq, before the variation of the pulse 
length is addressed. A global analysis of the dy- 
namics will be given in terms of the following quan- 
tities: (i) the energy absorbed by the molecule at 
the end of the pulse: AE^oi = (*(ip)|-ffmoi|*(ip))- 
(^'o|i?moi|^o)i (ii) the maximum energy change of the 
bath modes reached during the considered time inter- 

val of 5 ps: AEi^::^^ ^ma^i{^'it)\H^\^ij{t))-{^o\H^\^o)), 
(iii) the maximum expectation value of the Ei modes' 
coordinate operators during the 5 ps time interval: 
Q('"'^'^)=max((^'(t)|Q|*(i))), and (iv) the population 
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FIG. 2. Classical picture of vibrational wave packet circulation triggered by excitation of the two degenerate Ei modes by 
means of a circularly polarized laser field. Shown are the displacement vectors for the superposition of the modes along a unit 
circle in steps of 45°. In the central panel the two-dimensional potential V(Q, q = 0) is shown (contour lines in steps of 2000 
cm~^ starting at 2000 cm~^ ). 



of the ground state at the end of the pulse: Pq = 
l(*(tp)l*o)p. 

A. Dependence on the Field Strength 

A global analysis of the dynamics in dependence on 
the field strength and for a pulse duration of tp — 500 fs 
is given in Tab. |IV[ Overall, we observe a monotonous 
decrease of the ground state population Pq with increas- 



ing field amplitude. For the highest amplitude used, the 
ground state is completely depopulated. At the same 
time the absorbed energy AE^^i increases, indicating vi- 
brational ladder climbing. In fact for Eq = 1.5 niEh/eao 
the absorbed energy is compatible with an appreciable si- 
multaneous population of the third vibrationally excited 
state in each of the two Q modes. Analysis of the popula- 
tion shows that in fact many overtone and combination 
transitions are excited with the {vi = 2,1^2 = 2) state 
having the largest population. For the more moderate 
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TABLE IV. Analysis of the laser-driven wave packet dynam- 
ics for a 500 fs pulse which is resonant to the fundamental 
transition of the Ei modes (cf. Eq. (|6|). 



Eo 


ASmol 


A p(max) 


Q(max) 


Po 


(mEh/eao) 


(cm-1 ) 


(cm-1 ) 


(aoVa.m.u) 




0.1 


71 


6 


0.05 


0.97 


0.2 


277 


18 


0.09 


0.87 


0.3 


621 


62 


0.13 


0.73 


0.4 


1099 


164 


0.18 


0.57 


0.5 


1708 


341 


0.22 


0.42 


1.5 


12733 


4553 


0.50 


0.00 



excitation with Eq — 0.5 mEh/eag the populations of 
the Q modes' states is about 20% for (0,1) and (0,1), 
8% for (1, 1) and 4 % for (2, 0) and (0, 2). 

The anharmonic coupling to the bath modes will in- 
crease along with the absorbed energy, i.e. if the wave 
packet explores configurations away from equilibrium. 
For moderate excitation conditions the extent of anhar- 
monic coupling is rather small, e.g. if the amount of 
absorbed energy corresponds to about one quantum of 
excitation in a CO mode {Eq = 0.5 mEh/eao) we find 

^£'moi/^-£'battr'' ~ effect of the bath becomes 

appreciable only under high excitation conditions, e.g.. 



(max) 
bath 



2.8. In 



for Eq = 1.5 mEh/eao we find AE^-ioi/ AE^^ 
Fig. [3]we show the expectation values of the Hamiltonian 
operators for the system modes, Hq, and the uncoupled 
bath modes, Hq, for two field strengths. As expected 
the energy of the Q modes rises during laser pulse ex- 
citation (ip = 500 fs) as seen in panel (a). For the 
lower field strength the expectation value subsequently 
oscillates around the attained value signaling reversible 
energy exchange with the bath modes. Indeed the corre- 
sponding oscillations in the bath mode energies are seen 
in panel (b) . Increasing the field strength, the energy ex- 
change with the bath is more pronounced and instead of 
regular oscillations the onset of equilibration among the 
bath modes is observed in panel (c). Still there is no ap- 
preciable coupling to lower frequency modes that could 
lead to an energy flow out of the selected Q modes. We 
notice from panels (b) and (c) that the low-frequency 
mode qi, which was selected because it had the largest 
reorganization energy in the lower-frequency part of the 
spectrum is not appreciably excited. Finally, we com- 
ment on the initial partitioning of energy among the bath 
modes. Inspecting Tab. Ill one might have anticipated 
that mode will dominate the bath dynamics. However, 
the fundamental frequency of the laser-excited Q modes 
is 2014 cm~^ , which is resonant to q2 — qi, but rather 
off- resonant with respect to q^. Hence, energy exchange 
between the Q and the q^ modes is not very efficient. 

The focus of the present investigation is on the wave 
packet circulation with respect to the Ei modes. Inspect- 
ing Tab. |IV| we notice that the maximum bond elonga- 
tion increases with the amount of absorbed energy. The 




1000 2000 3000 4000 5000 

t[fs] 

FIG. 3. (a) Expectation values of Hq for a 500 fs excitation 
with field strength Eo = 0.5 mEh/eao (lower curve) and 1.5 
mEh/eao (upper curve). The respective expectation values of 
-ffq for the five bath modes are shown in panels (b) [Eo = 
0.5 mEh/eao) and (c) (JSq = 1.5 mEh/eao). With increasing 
energy at about 1000 fs the curves correspond to q\, qz,, qs, 
52, and 54 in panel (c), whereas in panel (b) §2 and §4 are 
interchanged (see also key in panel (c)). 



value that is reached can be compared, e.g., with the vari- 
ance of the groimd state coordinate distribution which is 
0.17 aoVa.m.u for the Q modes. Thus for the largest field 
strength we observe that the bond elongation is about 
three times this ground state variance. A more detailed 
picture is provided by inspecting the coordinate oper- 
ator expectation values, Qi{t) = {^{t)\Qi\^{t)), which 
are shown in Fig. [4] In the parametric plot the expecta- 
tion values trace a clockwise circulating trajectory on the 
potential V^(Q, q = 0). For small field strengths the am- 
plitude of this circulation increases during the excitation 
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TABLE V. Analysis of the laser-driven wave packet dynamics 
for pulses of different duration and amplitude. The integrated 
field envelope is identical in all cases. 





Eo 




A p(max) 


Q{max) 


Po 


(fs) 


(mEh/eao) 


(cm-^ ) 


(cm-1 ) 


(aoVa.m.u) 




100 


2.5 


1812 


374 


0.23 


0.40 


500 


0.5 


1708 


341 


0.22 


0.42 


1000 


0.25 


1441 


273 


0.20 


0.47 


1500 


0.17 


1146 


205 


0.18 


0.54 



period of 500 fs (upper and middle panel of Fig. |4]) . For 
the strongest field used (lower panel of Fig. |4]) the ampli- 
tude of the circulation starts to decrease already during 
the interaction with the pulse. Furthermore, after the 
excitation pulse is over the behavior shows a strong de- 
pendence on the field amplitude (right column of Fig. |4]) . 
Judging from the amount of energy redistribution given 
in Tab. IV the decrease of Qi{t) can only in part be an 
effect of the anharmonic coupling to the bath. Instead it 
is the wave packet dispersion that is mostly responsible 
for this decrease of Qi{t). This is shown in Fig. [s] where 
snapshots of the reduced density are plotted for the case 
of Eq = 0.5 mEh/eao- One notices that the wave packet 
stays compact in the excitation period only, but then 
starts to disperse. Eventually, the wave packet is rather 
delocalized even though it is still circulating. However, 
due to the delocalization negative and positive coordinate 
values contribute to the expectation values along the Q 
axes and therefore the expectation values become smaller 
and cannot refiect the wave packet circulation anymore. 
With increasing excitation more eigenstates contribute to 
the wave packet such that the dispersion effect becomes 
more pronounced and sets in already during the pulse. 



B. Dependence on the Pulse Duration 




t(ps) 



Qi 0.05^ 

OX 
-0.05, ' 



3.0 

2.5 1 
2.0[ 



:^05 1^ 



1.5j. 
1.0; 




t(ps) 




t(ps) 




FIG. 4. Parametric plots of (Qi, Q2) (in aoVa.m.u) for a 500 
fs laser pulse excitation and a field strength Eq (in mEh/eao) 
of 0.2 (upper row), 0.5 (middle row) and 1.5 (bottom row). 
(Notice the different axes scales). 



In the following we investigate the dependence of the 
laser-driven wave packet dynamics on the duration of the 
excitation pulse. As a reference case we have used t-p — 
500 fs and Eq = 0.5 mEh/eag. For a given duration tp 
the field amplitudes have been chosen such as to keep 
the integrated field envelope at the same value as in the 
reference case. Note that we are not dealing with a two 
level system and that the pulses are no 7r-pulses. The 
pulse parameters are summarized in Tab. [V| 

Table |V] also contains the various expectation values, 
which give a global characterization of the dynamics. 
First, we notice that the different pulse durations give 
comparable depopulations of the vibrational ground state 
between 40 and 54%. The amount of energy absorbed by 
the molecule decreases with the pulse duration and is 
largest for the shortest pulse {tp ~ 100 fs), although the 
difference to the reference pulse is only 6%. At the same 
time the excitation of bath modes decreases from a con- 



tribution to the total energy of 21% for tp — 100 fs down 
to 18% for tp = 1500. Moreover, the maximum radius 
of circulation of the expectation value of the Q coordi- 
nates is reduced slightly. Inspecting the time dependence 
of the Qi (t) in Fig. [6] we notice that again we have an 
initial period, where the wave packet stays rather com- 
pact such that the coordinate expectation values give a 
reasonable picture of the moving wave packet. This is fol- 
lowed by a period where the wave packet disperses such 
that the dynamics cannot be understood from the expec- 
tation values. As far as the pulse duration dependence 
is concerned the most striking observation from Fig. |6]is 
that the time span of compact wave packet propagation 
can be prolongated by increasing the parameter tp. In 
passing we note that this does also imply that the time 
span of circulation at about the full radius is longer, but 
of course not identical to tp. For tp = 500 fs and 1500 
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Os (ao(a.m.u)^'^) 



Oi (ao(a.m.u)''^) 




Qg (ao(a.m.u) 



1/2, 



-1 1 

Qi (ao(a.m.u)^'2) 



t=3ps 




A-1 

Qa (ao(a.m.u)^'^) 



Oi (ao(a.m.u)'''') 



FIG. 5. Snapshots of the time evolution of the reduced density 
with respect to the Q modes for laser-driving with a pulse 
having the parameters tp = 500 fs, Eq = 0.5 mEh/eao, and 
uj/2ttc = 2014 cm"^ . 



fs the expectation values Qi{t) are within 90% of their 
maximum value during the interval [450 : 800] fs and 
[950 : 1650] fs, respectively. Since the anharmonic cou- 
pling is still small in the considered range of the potential 
energy surface corresponding to about one to two quanta 
of excitation of the Ei modes, intramolecular energy re- 
distribution doesn't destroy wave packet circulation even 
for the longest pulse duration of 1.5 ps and subsequent 
3.5 ps free evolution. 



IV. SUMMARY 

We have demonstrated laser-driven wave packet 
circulation with respect to the two degenerate Ei 
carbonyl stretching vibrations in a pre-oriented 
Mn2(CO)io model. The wave packet circulation 
corresponds to a vibrational excitation moving around 



t(fs) 



t(ps) 




0.2-0.2 O2 



t(ps) 




XI ' — ^0-1 



t(ps) 




t(ps) 




Q, O^^S^ 

X\ ^^^-_0.1 



t(ps) 




FIG. 6. Parametric plots of (Qi,Q2) (in ao\/a.m.u) for dif- 
ferent pulse durations. Upper row: tp = 100 fs, Eo = 2.5 
mEh/eao, middle row: tp = 1000 fs, Eq = 0.25 mEh/eao, 
bottom row: tp — 1500 fs, Eo — 0.17 mEh/eao. (Notice the 
different axes scales). 



the symmetry axis of the system. It has been shown 
that a classical picture holds approximately during the 
interaction with the excitation laser pulse. Later on, 
appreciable wave packet dispersion sets in such that 
the classical picture becomes meaningless, although the 
wave packet still circulates around the minimum of the 
potential energy surface. 

The possibility of wave packet circulation upon exci- 
tation of degenerate vibrational modes with circularly 
polarized laser light could have been anticipated fr om 
the work of J. Manz et al. on triatomic molecules.l^ 
Therefore, the most important result of the present study 
concerns the robustness of wave packet circulation with 
respect to intramolecular vibrational energy redistribu- 
tion in a polyatomic molecule. To proof this we have 
considered the coupling of the two Ei mode to other in- 
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tramolecular degrees of freedom, which have been de- 
scribed in harmonic approximation. This resulted in a 
seven-dimensional model for which the time-dependent 
Schrodinger equation has been solved using the MCTDH 
approach. For moderate excitation conditions of one or 
two quanta in the Ei modes, at most about 20 % of the 
absorbed energy is disposed into the bath modes within 
a time span of 5 picoseconds. This finding appears to be 
encouraging and should stimulate experimental verifica- 
tion. 
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